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A theoretical description of the sequence of magnetic phases in CosTeOe is presented. The 
strongly first-order character of the transition to the commensurate multiferroic ground state, in- 
duced by coupled order parameters corresponding to different wavevectors, is related to a large 
magnetoelastic effect with an exchange energy critically sensitive to the interatomic spacing. The 
monoclinic magnetic symmetry C2' of the multiferroic phase permits spontaneous polarization and 
magnetization as well as the linear magnetoelectric effect. The existence of weakly ferromagnetic 
domains is verified experimentally by second harmonic generation measurements. 
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At variance with structural transitions which alter the 
lengths and orientations of the chemical bonds, the spin 
ordering occurring in magnetic phases has in most cases 
a negligible effect on the structural lattice. An im- 
portant exception is represented by the class of mag- 
netostructural transitions occurring in multiferroic com- 
pounds in which the magnetic ordering in the multifer- 
roic phase induces simultaneously a change in the atomic 
structure, w hich permits the emergence of a spontaneous 
polarization.ll^l However, the measured changes of lat- 
tice parameters found in the multiferroic phases are gen- 
erally small, i.e. of the order of 10~^ A, and do not af- 
fect the second-order character of the transitions to these 
phases .'^^ Here, we describe theoretically the sequence 
of phases recently reported in CoaTeOj^EJ jj^ which a 
strongly first order transition, characterized by substan- 
tial discontinuities of the lattice parameters and a re- 
markable delta-shape peak of the specific heat, yields 
a multiferroic ground state displaying magnetoelectric 
properties. The observed structural changes are related 
to the coupling between the magnetic order-parameters 
involved at the transition, which correspond to different 
propagation wave-vectors, in contrast with the standard 
situation found in multiferroic transitions where the cou- 
pled order-parameters generally pertain to the same k- 
vector.!^ 



Neutron powder diffraction studies^ show that below 
the paramagnetic phase described by the space group 
Gp = C2/cl' the CosTeOe undergoes a sequence of 
three antiferromagnetic phases, summarized in Fig. [Ipl 
They are associated with three different k- vectors of the 
centred monoclinic Brillouin-zone: ki = (0,0.480,0.055), 
h = (0,0,0) and fca = (0,1/2,1/4). 

L The incommensurate Phase I, which emerges at 
Tjv — 26 K, shows the coexistence of fci and fc2- The tran- 
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FIG. 1. Sequence of phases observed in CosTeOe from Ref. |B] 

sition order-parameter transforming as the irreducible 
representation (IR) Ti{ki) (see Table hP^ has four com- 



ponents (?7i = pie'^'S rjl = Pie '"^ 772 = P2e'''^ 772 = 
P2e~ ) which yield the free-energy 

Fi{pi,P2) = ^{pl + pl) + ^{pt + pt) + ^plpl (1) 

Minimizing Fi leads to two possibly stable states having 
the same point-group symmetry 2/ml' as the paramag- 
netic phase, and corresponding to the equilibrium condi- 
tions (pi ^ 0, P2 — 0) and (pi = ±/02) ■ The point group 
symmetry 2/ml' of phase I and its incommensurate char- 
acter are preserved when considering a coupling with the 
order-parameter associated with fc2- 

H. The commensurate phase H appearing at Ti = 
21.1 K below a second-order transition corresponds to 
the single wave- vector fc2- It remains stable in a nar- 
row interval of temperature down to T2 — 17.4 K. The 
one-dimensional IRs Fi — r4 at the centre of the mono- 
clinic Brillouin zonc"'^^ induce respectively the magnetic 



TABLE I. Generators of the active irreducible representations of the paramagnetic space-group C2/cl' associated with the 
wave- vectors fci, fcg and fca in CosTeOe. Columns matrices represent diagonal matrices. T is the time-reversal operator, 
e = exp (ifcf c), uj = exp (\k\b/2) 
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symmetries C2/c (Fi), C27c' (Fa), C2/c' (Fg) and C27c 
(r4). The magnetic structure proposed by Ivanov et 
al.l^ from neutron data coincides with the C2'/c mag- 
netic group which is therefore associated with a one- 
dimensional order-parameter, denoted C, hereafter, corre- 
sponding to the single equilibrium state which minimizes 
below Ti the canonical free-energy: 



^2(C) = fc^ + ^C^ 



(2) 



In absence of applied fields the magnetic symmetry of 
the phase does not allow the emergence of spontaneous 
polarization or magnetization components. 

III. At T2 = 17.4 K a commensurate phase III arises in 
which the neutron diffraction pattern corresponding to ^2 
persists coexisting with magnetic peaks associated with 
the commensurate wave- vector k^. Both sets of reflec- 
tions are observed in the whole range of stability of the 
phase down to 1.6 K.^ k^ is in general a position inside 
the Brillouin-zone, corresponding to a four-dimensional 
IR of Gp, denoted Ti(fc3), whose matrices are listed in Ta- 
ble |T] Keeping for the four order-parameter components 
the same notation {-qi = pie^^^\i = 1,2) as for phase I 
the transition free-energy reads 



Fz{pi,p2, 01,62) = Fi{pi,p2) 



P3 



(Pi COS 46*1 + P2 cos 46*2 ) 



(3) 



^p?p^sin2(0i+02) + 



which differs from Fi by the /^a and (3^ lock-in invari- 
ants. Phase III results from the coupling of the order- 
parameters rii{k^) and C(^2) corresponding to the total 
free-energy 

FT = F:,{p,,e,) + F2{0-\e{pl+pl) , (4) 

where the 5-term represents the lowest-degree coupling 
between the two order parameters. Table In] lists the 



symmetries and equilibrium conditions of the seven pos- 
sibly stable magnetic phases resulting from the minimiza- 
tion of Ft- The phase with magnetic symmetry C2' and 
a sixteen-fold multiplication of the primitive paramag- 
netic unit cell, shown in Fig. ^a), coincides unambigu- 
ously with the reported neutron diffraction observationsSl 
in phase III. It corresponds to the equilibrium values of 
the order-parameters C 7^ ^j Pi = P2 = Pe,di = O2 = 0^ 
which allow emergence of a spontaneous polarization Py 

and a spontaneous weak magnetization M = {Mx,Mz), 
with, respectively, two ferroelectric and two weakly fer- 
romagnetic domains (Fig. [2Fb)). 



b) 






FIG. 2. (a) Paramagnetic unit-cell embedded into the 16-fold 
unit-cell of the multiferroic phase III of CosTeOe. (b) Ferro- 
electric and weak ferromagnetic domains in phase III. In the 
text we use Cartesian instead of the monoclinic coordinates 
according to a; ~ a, j/ ~ 6, z ~ c. 



TABLE II. Seven possible choices (a) of the magnetic space groups (b) derived from the minimization of Ft (C,Pi, Si) in Eq. (Ml. 
(c) Equilibrium values of the order parameters, (d) Basic translations of the conventional monoclinic or triclinic unit cells, (e) 
Multiplicity of the volume of the primitive paramagnetic unit cell, (f) Origin of the coordinates. 



(a) 
1 



(b) 
PI' 

PI' 



3 


PI 


4 


C2' 


5 


PI' 


6 


PI' 


7 


PI 



(c) 

C/0,pi /0,p2 = 0,6»i =0 

C/0,pi /0,p2 = 0,6»i =7r/4 

C/0,pi/0,p2 = 
C 7^0, pi =p2,ei =02 
C7^0,pi 7^ 0, p2 / 0, 6)i = 0, 02 



-7r/4 



C/0,pi /0,p2/0,6)i=7r/2,6'2 = -7r/4 



C/0,pi/0,p2/0 



(d) 



1 




(c) 


(f) 






4 


(0,1/2,0) 


(-1,0,0) 








(1/2,3/2,1) 




4 


(1/4,3/4,0) 


(1/2,1/2,-1) 












4 


(0,0,0) 


0,4), (0,2,0), (- 


-1,0,0) 


16 


(0,5/8,5/4) 






8 


(0,1/2,0) 


(0,2,0) 








(1,0,0) 




8 


(1/4,1/4,-1/2) 


(0,1,-2) 












8 


(0,0,0) 



Figure l3] shows the distribution of magnetic domains 
in the monoclinic xz plane. The image was gained by 
optical second harmonic generation (SHG) as described 
in Ref. 5. The orientation of the domain walls along 
an arbitrary direction in the xz plane further confirms 
the 2' symmetry. This symmetry is also consistent with 
the presence of the Xxxx and Xzzz components of the 
SHG susceptibility tensor.El Note, however, that it dif- 
fers from the symmetry m assumed in Ref. El In Ref. [5] 
only i-tensor components were considered as origin of the 
SHG signal since in magnetically induced ferroelectrics 
like MnW04 the SHG signal is always linearly coupling 
to the spontaneous polarization. In contrast, the SHG 
signal leading to Fig. |3] revealed that SHG in GTO is 
related to c-tensor components reproducing the weakly 
ferromagnetic order. SHG with Xxxx 7^ 0, Xyyy = 0, and 
Xzzz 7^ as c-type susceptibilities leads to the magnetic 
symmetry 2'. Unfortunately, a recent discussion of our 
SHG data in Ref. [TT] is still based on the assumption of 
SHG coupling to the electric polarization which lead to 
results inconsistent with the ones reported in this work. 

Since the transitions to phases I and III result from 
the coupling of two order-parameters corresponding to 
distinct fc-vectors they display necessarily a first-order 
charactepi2J following the triggering mechanism proposed 
by HolakovskjU^ in which one order-parameter triggers 
the onset of another order-parameter across the first- 
order discontinuity. The mechanism requires taking into 
account a sixth-degree invariant of the "triggered" order- 
parameter (^(pi + P2)^) in the total free-energy Ft- Un- 



der the conditions S > 



|ffi+fe|Ai 
2 



1/2 



and Ai - 



iS^ 



< 



l^i+fcl 

^ y ^ J , the phase corresponding to the sole "trig- 
gering" order-parameter ^ becomes unstable with re- 
spect to a phase in which both order parameters, C o,nd 
rj ~ (pi,9i), are non-zero. In this phase C is frozen and 
rj determines the symmetry breaking process. The tran- 
sition can be shown to occur discontinuously at a higher 




FIG. 3. SHG image of a polished single-crystal platelet 
(thickness 50pm) of CosTeOe at T = 5K. Diagonal paral- 
lel stripes correspond to magnetic domain walls, whereas the 
gradual ringlike oscillation of brightness is caused by interfer- 
ence of the laser light in the sample. The incident light at 
E = 1.29 eV was propagating parallel to the crystallographic 
y axis and polarized parallel to the direction of the domain 
walls, while the SHG light was polarized perpendicular to the 
domain walls. 



temperature than the transition temperature at which 
a phase with ( — 0,pi ^ would appear^ Accord- 
ingly, the triggering process, which is activated in the 
region of phase coexistence preceding the transition at 
T2, is due to a large negative value of the interaction 
term for the coupling between the two order-parameters, 
which determines the value of the coupling coefficient S. 
Note that the first-order character of the transition to 
the multiferroic phase at T2 is confirmed by the strong 
discontinuities (of the order of 10~^ A) observed in the 
lattice paramete rs^ an d by a remarkably sharp peak of 
the specific heat.l^EES Jn contrast, the transition from 
the paramagnetic to the antiferromagnetic phase I is 



weakly first-order with almost negligible lattice discon- 
tinuities and a standard specific heat anomaly!^ The 
transition from phase I to phase II which involves a sin- 
gle order-parameter has typical second-order transition 
features with no noticeable discontinuity of the lattice 
parameterspEl 

The dielectric contribution to the free-energy F^, = 

p2 

—vPyC^iplpl cos 2{9i + 62)) + 03- yields the equihbrium 

yy 

value of Py below T2 : 



Py = iye'yX' pi cos A9, 



(5) 



At T — T2, Py undergoes an upward discontinuity, im- 
posed by the first-order character of the transition. On 
further cooling it increases as sa (T2 — T)^, since the C 
order-parameter is frozen in phase III. A similar temper- 
ature dependence holds for the spontaneous magnetiza- 
tion components M^ and M^. From the spontaneous 
magnetic contribution to the free-energy in phase III 



Fm = PuMuC,plplcoa2{6i 
one gets 



O2 






(with u = x,z), 



M„, 



-xL^uCPe cos 46*6 



(6) 



as the equilibrium value below T2. Application of mag- 
netic fields along x or z yields a renormalization of the 
transition temperature according to 



T2{H,,)=T2{Q)~xluPua^\''Hl 



(7) 



Where a — aQ{T ~ 72(0)) is the coefficient of p^ in 
Ft- For pu > the transition temperature is low- 
ered under application of a magnetic field and T2(0) — 
T2{Hu) increases quadratically with H^, as observed in 
CoaTeOg under H^ field. The magnetic susceptibility 
components Xuu = jj^ are obtained by minimizing the 



field-induced contribution to the free-energy Fm{Hu) 
""^"^ ' 2-rr HuM-a- This reveals 



p^m^h^CpI 



Xv. 



X°u(l-MuC^Pe) 



(8) 



Therefore the discontinuous jump of the order- 
parameter Pe(T2) coincides with a drop of XuuiT2) the 
magnitude of which decreases with increasing field up 
to a threshold field corresponding to Pe{H^) — ^ j — 

above which Xuu undergoes an upward discontinuity at 
T2iHu). This behaviour is verified experimentalljffl for 
H^ with a threshold field H^ « 12 T. The slight de- 
crease with temperature observed for Xzz below T2 in- 
dicates a slight increase of the order-parameter within 
the multiferroic phase, with an almost step-like depen- 
dence on temperature across the transition. As a conse- 
quence of this Heaviside-like behaviour the specific heat 
C which is proportional to the order-parameter deriv ative 
displays a delta-shape like behaviour across T2 ! ^ * ^ * ^^ ^ This 
remarkable property of the first-order multiferroic tran- 
sition in CoaTeOg, which has been previously obs erved 
at the first-order ferromagnetic transition in Fe2P,'^^'^ 



is reminiscent of structural transitions having a recon- 
structive mechanisnP^ as observed, for example, at the 
fcc-hcp transition in cobalt .1^ However, in CosTeOg the 
symmetry-breaking mechanism involved at T2 is not of 
the reconstructive type since the C2' symmetry of phase 
III is group-subgroup related to the C2'/c symmetry of 
phase II. Therefore the strong discontinuities reported 
for the lattice parameters^ at T2 should correspond to 
a strong magnetoelastic coupling with an exchange en- 
ergy critically sensitive to the interatomic spacing. This 
is in agreement with the large spin lattice coupling de- 
duced by Her et al.l^ from the magnetic hysteresis curves 
and with the direct exchange pathway corresponding to 
the shorter Co-Co distances existing in phase III (see 
Fig. 12 in Ref. [5]) for the neighbouring atoms Co(5)- 
Co(5), Co(2)-Co(5), Co(3)-Co(3) and Co(4)-Co(4). The 
strong magnetoelastic effect is favored by the coupling 
between the order-parameters C and rji. Here, the already 
existing antiferromagnetic order-parameter C triggers the 
emergence of the order parameter rji which is reflected by 
the transition discontinuity. 

The respective monoclinic symmetries of phases II and 
III permit a variety of magnetoelectric effects under ap- 
plied magnetic or electric fields. For instance, in phase II 
applying a magnetic field Hy induces contributions by Px 
and Pz to the polarization which vary as Px^z ~ C^y In 
phase HI one has Px^z ~ (,ptHyi i-^-, at constant temper- 
ature Px and Pz increase linearly with Hy in both phases, 
while at constant field they increase with temperature as 
(Ti - T)i/2 in phase I and as {T2 - T) in phase HI. Con- 
versely, applying the magnetic field Hx or Hz in phase 
II induces a polarization polarization Py « C,Hx,z- In 
phase HI an additional contribution to the spontaneous 
polarization Py, namely APj, « C,p\Hx^z is generated. 
Reversed magnetoelectric effects should also be observed 
in phases II and HI with the onset of an induced magneti- 
zation My depending linearly on electric fields Ex or Ez ■ 
Additional contributions to the Mx and Mz components 
of the spontaneous magnetization are generated under 
an electric field Ey . Accordingly, our theoretical analysis 
suggests that magnetic-field-induced polarization compo- 
nents Px{Hx), Px{Hz) and Pz{Hx) are not allowed. Their 
observation was reported earliep^and should be related to 
a misorientation of the sample admixing a y component 
to the X axis. The possibility for such a misorientation 
was indeed pointed out by the authors of Ref. ,5, 

In summary, the observed properties of the sequence 
of three phases reported in CoaTeOg have been described 
theoretically. The most striking feature of this phase se- 
quence is the coexistence of propagation vectors in the 
incommensurate and commensurate multiferroic phases, 
the Brillouin zone-centre ^2 vector, which persists in all 
three phases, coupling successively with the incommen- 
surate fci-vector in phase I, and with its lock-in commen- 
surate variant fcs in phase HI. Another remarkable prop- 
erty is the strongly first-order-character of the multifer- 
roic transition which relates both to the triggering mech- 
anism coupling the ^2 and k^ related order-parameters. 



and to a strong magnetoelastic coupling. The sharp peak 
of the specific heat is shown to be consistent with the al- 
most constant value of the spontaneous magnetization in 
the multiferroic state. Furthermore, a number of experi- 
mental observations of Ref. i5]have been scrutinized: The 
monoclinic symmetry ofphase III is 2' (instead of the 
m as previously proposetpl) . This permits a spontaneous 
weak magnetization {M^, M^), which is found to be in 
agreement with domain structures observed in SHG mea- 
surements. It also allows a spontaneous polarization Py 
which has not been investigated before whereas contri- 
butions Px^z 7^ are no longer expected. The magneto- 
electric effects that should exist in CoaTeOe have been 
worked out theoretically. A verification of these predic- 
tions, supported by the application of both magnetic and 
electric fields, is necessary for confirming the validity of 
the theoretical description presented in this article. 

At last we emphasize that our theoretical analysis and 
conclusions differ in an essential manner from the sym- 
metry analysis proposed in Ref. [TT] for two reasons, (i) 
Scrutinization of the data reported in Ref. [5| lead to a 
revision of the magnetic symmetry and the direction of 



electric polarization. Since this is a very recent result 
it could be taken into account in the present work but 
not in Ref. [TTJ (ii) In Ref. [TT] mainly one dimensional 

order-parameters associated with the wave-vector ^2 are 
considered. In contrast, we analyze the symmetries and 
physical properties of phases I and III as resulting from 
the coupling of order-parameters associated respectively 
with the wave-vectors (fci, k2) and (fc2, k^), in agree- 
ment with the neutron diffraction data reported in Ref. [HI 
In particular the field- induced component Pz{Hx) dis- 
cussed in Ref. II H is shown to be absent in our description, 
whereas we predict the existence of a single spontaneous 
(zero- field) polarization component Py. 
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